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ABSTRACT 


Noise  power  on  the  short  period  pendulum  and  strain 
instruments  has  been  compared  at  WMO  and  HNME.  In  the  re¬ 
sponse  band  of  the  short  period  instru?ients ,  .here  are 
approximately  four  decades  of  relative  power  cl  HNME,  and 
two  decades  at  WMO.  In  one  case  analyzed  the  noise  power 
at  HNME  is  about  2%  times  the  noise  power  at  WMO  at  1  Hz. 

Coherence  was  estimated  as  a  measure  of  the  existence 
of  a  linear  transfer  function  between  the  vertical  strain 
and  the  vertical  pendulum.  At  HNME  there  exists  a  well- 
defined  linear  transfer  function  in  the  band  0.1  to  0.6  Hz. 
The  low  coherence  at  higher  frequencies  at  HNME  and  through¬ 
out  the  0.1  to  3.0  Hz  band  at  WMO  rule  out  the  possibility 
of  a  linear  transfer  function  between  the  vertical  strain 
and  vertical  pendulum. 

Multiple  coherence,  and  rotation  of  the  horizontal 
short  period  seismograms  were  used  to  infer  the  existence  of 
unidirectional  noise  components.  At  HNME  about  80  percent 
of  the  noise  power  is  unidirectional  in  the  0.3  to  0.4  Hz 
band.  At  WMO  the  noise  field  has  no  apparent  single  pre¬ 
ferred  direction.  Preliminary  work  on  rotation  of  the 
horizontal  strain  instruments  at  WMO  is  presented.  The  tech¬ 
nique  uses  a  third  horizontal  instrument  to  establish  the 
shear  component  of  strain. 

Examples  of  detailed  analysis  to  establish  wavetype  using 
the  vertical  strain,  together  with  the  vertical  pendulum  and 
horizontal  pendulum  instruments  are  included. 
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INTRODUCTION 


Benioff  (1935)  described  a  linear  strain  seismograph 
and  summarized  early  work  on  earth  strain;  Benioff  and 
Gutenburg  (195?)  compared  the  response  of  strain  and  pen¬ 
dulum  seismographs  to  earthquake  generated  Rayleigh  waves. 
Romney  (1964)  and  Benioff  (1962)  from  somewhat  different 
aspects  considered  the  response  to  surface  waves  and  P-waves , 
of  combinations  of  pendulum  and  strain  seismographs.  Smith 
(1966)  applied  strain  seismographs  to  analysis  of  free 
oscillations  of  the  earth  excited  by  large  earthquakes. 
Shopland  (1966)  described  the  strain  installation  of  WMO. 

Press  (1965)  exhibited  two  examples  of  apparent  permanent 
changes  in  the  strain  associated  with  earthquake  events. 

Gupta  (1966)  published  a  theoretical  study  of  strain  to  pen¬ 
dulum  seismograph  respor^ a  ratios. 

The  above  cited  papers  taken  together  with  a  number  of 
reports  published  by  Teledyne  Industries,  Geotech  Division, 
constitute  the  background  for  this  study,  which  attempts  to 
describe  the  noise  field  at  several  sites  and  to  develop  meth¬ 
ods  of  signal  enhancement  using  strain  seismometers. 

Data  from  two  sites  are  examined  in  this  report.  These 
are  WMO  (Wichita  Mountains  Observatory)  and  HNME  (The  Long 
Range  Seismic  Measurement  station  at  Houlton,  Maine).  Both 
sites  have  long  and  short  period  pendulum  instruments;  HNME 
has  a  vertical  strain  instrument;  and  WMO  has  orthogonal  strain 
instruments.  The  data  used  in  preparation  of  this  report  is 
listed  in  Appendix  I. 
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TECHNICAL  T-xSCUSSION 


Site  geology 

WMO  is  located  in  the  Wichita  Mountains  of  West  Central 
Oklahoma.  The  Wichita  Mountains  are  the  outcropping  portion 
of  a  major  granitic  thrust  from  the  south.  The  Anadarko 
basin  north  of  the  Wichitas  is  asymmetric;  its  deepest  area 
being  near  the  mountains.  The  north  flank  of  the  Wichitas 
exhibits  evidence  of  overthrusting  and  complex  folding. 

HNME  is  located  in  east  central  Maine  near  the  New 
Brunswick  border.  The  outcropping  rocks  are  Silurian,  and 
consist  of  interbedded  marine  sediments,  tuffs,  and  lavas. 

The  general  region  is  characterized  by  complex  intrusions 
and  folding.  The  youngest  rocks  in  the  region  are  probably 
Paleozoic.  It  is  approximately  150  km  from  the  Atlantic  coast 
line. 

Noise  power  compa^'son 

Figures  1  and  2  compare  relative  noise  power  spectra  as 
recorded  on  the  vertical  short  period  pendulum  instruments  at 
WMO  and  HNME.  Calibration  of  Figures  1  and  2  shows  the  noise 

rs 

level  to  be  0.52  (my/sec)  /Hz  at  1.0  Hz  at  WMO  and  1.29  (my/sec)2/ 
Hz  at  1.0  Hz  at  HNME.  The  noise  at  WMO  includes  about  two  orders 
of  magnitude  of  relative  power;  at  HNME  it  covers  about  four 
orders  of  magnitude.  These  plots  have  not  been  adjusted  for 
system  response.  Both  pendulum  instruments  show  a  narrow  peak  in 
the  noise  at  higher  frequencies,  2.6  cps  at  HNME  and  2.0  cps  at 
WMO.  Figure  3  is  noise  power  at  HNME,  with  less  smoothing, 
included  to  illustrate  better  the  noise  peak  at  2.6  cps.  These 
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plots  of  short  period  power  illustrate  the  noise  field  in 
which  the  strain  instruments  have  been  installed. 

Figures  4  and  5  are  comparable  plots  of  seismic  noise 
as  recorded  by  the  vertical  strain  and  vertical  pendulum  at 
WMO.  The  higher  relative  power  of  the  2.0  cps  noise  on  the 
pendulum  as  compared  with  the  strain  instrument  is  typical 
of  a  number  of  samples  observed.  This  is  consistent  with  the 
hypothesis  that  the  noise  is  nearly  vertical  P  waves  which 
are  not  picked  up  on  vertical  strain  instrument.  This  sample  also 
contains  a  peak  in  the  noise  at  0.95  cps  that  has  higher 
relative  power  on  the  pendulum  instrument. 

Approximate  system  noise  level  on  the  power  spectra 
plots  has  been  determined  by  observing  the  level  at  which 
the  noise  is  white  in  the  5.0  to  10.0  Hz  region;  the  assump¬ 
tion  being  that  system  noise  is  white  throughout  the  range 
of  frequencies  analyzed. 

Figures  6  and  7  compare  the  vertical  pendulum  and  verti¬ 
cal  strain  at  HNME.  These  plots  show  the  increase  in  per¬ 
centage  variation  of  power  from  0.4  Hz  to  3.5  Hz  and  also  show 
the  greater  relative  response  of  the  pendulum  instrument  to 
the  2.6  Hz  noise.  The  energy  at  this  frequency  will  be 
examined  in  some  detail  in  a  later  section  of  this  report. 

Squared  coherence  comparison 

Figures  8  and  9  are  plots  of  th.e  squared  coherence  between 
the  vertical  strain  and  vertical  pendulum  a*-  HNME  and  WMO. 

These  are  typical  of  a  number  of  samples  of  varying  length  and 
varying  amounts  of  smoothing  that  were  obtained. 
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The  squared  coherence  is.  defined  by 


<2 

12 


(f) 


cf ) 


in  which  ot^2  is  the  sum  of  the  squares  of  the  cospectrum 
and  the  quadspectrum  between  the  two  noise  processes,  and 

and  r22(f )  are  the  power  spectra  of  the  two  processes, 
see  Jenkins  and  Watts  (1968)  p.  352. 


At  HNME  the  squared  coherence  between  the  vertical  strain 
and  vertical  pendulum  is  greater  than  0.9  between  0.15  and  0.4  Hz 
and  greater  than  0.7  out  t n  0.6  Hz.  The  HNME  plot  (Figure  8)  is 
interpreted  to  mean  that  there  exists  a  well-defined  linear  trans 
fer  function  between  the  vertical  strain,  and  vertical  pendulum 
in  the  0.15  to  0.6  Hz  band;  the  decrease  in  squared  coherence  for 
frequencies  >0.6  Hz  can  be  caused  either  by  added  random  noise 
on  one  of  the  instruments  or  by  a  mixture  of  two  or  more  linear 
processes . 


Each  of  the  plots  of  coherence,  and  later  of  multiple  co¬ 
herence,  includes  a  horizontal  dashed  line,  which  marks  the  95 
percent  confidence  limit  that  the  true  squared  coherence  is  zero. 
This  limit  is  obtained  by  assuming  two  independent  time  series 
and  noting  that  the  estimated  squared  coherence  approximately 
follows  a  Fisher  F-dis tribution.  Substitution  of  the  number  of 
degrees  of  freedom  and  use  of  tabulated  data  produces  a  quadratic 
which  may  be  solved  for  the  95  percent  confidence  limit.  Experi¬ 
ence  using  these  results  on  a  number  of  samples  of  varying 
lengths  with  a  variety  of  degrees  of  freedom,  indicates  the 
validity  of  the  95  percent  confidence  limit.  This  can  be  seen, 
for  example,  in  Figure  8.  (Jenkins  and  Watts  (.1968  )  p.  433). 

It  is  also  possible  to  .compute  confidence  limits  of  the 


observed  value  of  the  squared  coherence.  The  essential 
elements  of  this  procedure  are  the  transformation  of  the 
coherence  using  a  Fisher  z-transformation. 


Y12(f)  =  arctanh  I|K12|] 

Y12  1S  approximately  normally  distributed.  (Jenkins  and  Watts, 
1968,  p.  379-380;  Cramer,  1946,  p.  398-401.)  Bendat  and 
Piersoll  (1966,  p.  214)  use  a  similar  formula  to  establish  an 
estimated  value  of  squared  coherence.  For  K^2  =  0.9  and  32 
degrees  of  freedom,  the  95  percent  confidence  limits  are 

0.83  <  K^2  <  0.94 

which  fits  the  observed  variation  of  Figure  3. 

Figure  9  is  a  typical  plot  of  the  observed  squared  co¬ 
herence  at  WMO .  The  95  percent  confidence  limit  is  marked  on 
this  plot  for  0.155  and  0.31  Hz.  There  can  be  no  well-defined 
linear  transformation  between  the  vertical  strain  and  vertical 
pendulum  at  WMO.  Inspection  of  the  raw  data,  however,  reveal 
short  segments  of  5.0  to  10.0  seconds  duration  during  which  a 
linear  transformation  exists.  A  sample  of  this  sort  will  be 
examined  in  detail  in  a  following  section. 

Estimation  of  the  transfer  function 

In  order  to  obtain  an  estimate  of  the  linear  relation 

between  the  vertical  strain  and  vertical  pendulum  a  model  of 
the  form 

X2(t)  =  |  h(u)  X±  (t-u)  du  +  Z (t ) 
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is  taken,  following  Jenkins  and  Watts  (1968,  p.  351),  if 
the  cross  covariance  of  the  output  is  taken  and  then  Fourie 
transformed 


ri,(f) 

H(f)  -  T77TF) 


11 


in  which  ri2(f)  1S  the  complex  cross  spectrum  of  the  output 
and  input  and  f^^(t)  is  the  input  autospectrum. 


Writing 


H(f )  =  G(f  )e+llf’(f  5  = 


A12(f)  -  i¥12(f) 


in  which 


G(f)  =  [An%(f)  +  yf0(f)]1/  2/rni  (f )  =  o10(f)/r„(f) 


12 


12 


11 


12 


11 


which  defines  a12(f),  and 


$(f)  =  arctan  - 


12 


T^rn 


The  squared  coherency  is  defined  as 


(f ) 


‘12  "  r^Tfl  f22(f) 


H(f ) 


=  ..JWSl  1/2  e^(f) 


then 


In  comparing  HNME  and  WMO,  two  cases  occur,  figure  10 
is  plot  of  <j>(f)  fo^  several  examples  of  HNME  data.  For 
f  <_  1.0  Hz,  <£(f)  -  tt/2;  thus 

r  2  2  * 

H(f )  -K,  0 Cf  )I~— •]  (-i ) 

l  1± 

At  WMO  a  number  of  computations,  of  which  Figure  9  is 
an  example,  show  that  k^2  is  less  than  the  95  percent  confi¬ 
dence  limits  that  the  squared  coherence  is  zero.  That  is  at 
WMO 


H(f )  =  0 

In  the  model  assumed  this  means  that 


x2(t)  =  Z(t) 


that  is  the  vertical  pendulum  looks  like  noise  when  compared 
with  the  vertical  strain. 


Returning  to  HNME,  the  values  in 


H(f ) 


-jK12(f) 


are  given  in  Table  1,  at  0.31  Hz  for  a  sequence  of  512  point 
samples  1000  pts  (50  seconds)  apart. 


Thus  for  ar.y  given  moment  one  can  expect  any  fixed  linear 
transfer  function  to  fail  to  cancel  some  7  percent  of  the  noise 
amplitude  and  in  extreme  cases  as  much  as  10  percent  of  the 
noise  at  0.312  Hz.  It  is  assumed  that  the  vertical  strain  output 
is  phase  shifted  90°  and  then  subtracted  from  the  vertical 
pendulum.  This  is  another  way  of  looking  at  the  fact  that  the 
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coherence  is  not  100  percent,  i.e.,  some  of  the  power  on  the' 
yertical  pendulum  cannot  be  predicted  by  the  vertical  strain. 

Noise  source  direction 


Squared  Multiple  Coherence 


Two  techniques  are  used  to  analyze  the  noise  source 
direction  at  WMO  and  HNME .  They  are  squared  multiple  co¬ 
herence,  and  rotation  of  the  horizontal  instruments. 


The  squared  multiple  coherence  is  defined  as 


2 

«(q+l)123...q 


-  ^  F (q  +  1 ) (qt 1) (f } 1 

^(q+lMq+l)^  |?qq(f } 


A 

where  r  is  the  spectral  matrix  of  the  stochastic 
For  the  case  of  two  inputs  (q=2)  and  one  output  as 
this  report,  the  above  formulation  becomes 


process . 
used  in 
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where,  for  example,  T 12  is  the  cross-spectrum  between  channels 
1  and  2  and  F-^  is  the  autospectrum  of  channel  1.  The  squared 
multiple  coherence  measures  the  proportion  of  the  output  channel 
that  can  be  predicted  from  the  inputs  (.Jenkins  and  Watts,  1968  , 
p.  487-488).  As  applied  here,  the  horizontal  instruments  are 
used  as  inputs,  and  the  vertical  instrument  is  the  output.  If 
the  squared  multiple  coherence  is  high,  one  would  expect  the 
noise  to  be  unidirectional  and  therefore,  to  find  a  linear 
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transformation  hetween  horizontal  and  vertical  component  . 
of  the  noise.  Low  coherence  can  be  caused  either  by  multi¬ 
directional  noise  or  by  the  absence  of  a  time  invariant 
linear  transfer  function  between  the  horizontal  and  vertical 
components  of  the  noise. 

Figures  11  through  14  are  examples  of  the  squared  mul¬ 
tiple  coherence  as  obtained  at  WMO  and  HNME.  Typically,  the 
squared  multiple  coherence  has  values  near  0.8  at  0.3  Hz  at 
HNME;  at  WMO  the  value  is  consistently  below  or  near  the  95 
percent  confidence  limit  that  the  true  multiple  coherence  is 
different  from  zero.  The  high  squared  multiple  coherence  at 

2.6  Hz  at  HNME  appears  to  be  a  special  case  and  will  be  examined 
in  detail. 

Horizontal  pendulum  rotation 

If  there  exist  short  segments  (5  to  10  seconds)  of  the 
noise  that  are  unidirectional;  then  by  rotation  of  the  horizon¬ 
tal  pendulums  so  that  one  component  is  aligned  with  the  noise 
source  direction  and  the  other  is  perpendicular  to  it,  the 
maximum  amplitude  should  occur  on  one  component  and  the  minimum 
amplitude  on  the  component  at  90  degrees  to  the  azimuth. 

Figures  15  through  18  are  examples  of  the  horizontal  pen¬ 
dulums  at  HNME  rotated  in  15  degrees  steps  starting  from  their 
installed  azimuths.  The  sample  chosen  is  the  same  as  that 
used  to  obtain  the  squared  multiple  coherence  of  Figure  11. 
Clearly ,  if  the  noise  has  two  or  more  source  directions  at  the 
same  time,  this  technique  will  not  identify  the  direction. 

1 igure  15  is  the  original  rotated  seismograms.  Figures  16 
through  18  are  the  rotated  seismograms  filtered  as  marked. 

On  each  of  the  filtered  sets,  segments  of  noise  that  appear  to 


be  unidirectional  are  enclosed  in  a  box.  The  seismograms 
filtered  0.1  to  0.5  Hz  appear  to  have  a  predominance  of  noise 
from  3°  or  183°  (the  technique  does  not  allow  one  to  dis¬ 
tinguish  sense).  The  0.5  to  1.2  Hz  band  (Figure  16)  has  no 
apparent  source  direction,  yet  some  scattered  unidirectional 
segments  occur.  The  band  1.2  to  2.4  Hz  (Figure  17)  shows 
that  much  of  the  2.6  Hz  energy  has  an  azimuthal  direction  of 
168°  or  348°.  This  phenomenon  will  be  examined  in  detail. 

Figures  19  through  22  have  the  same  format  as  Figures  15 
through  18;  the  data  is  now  from  WMO.  The  same  technique  has 
been  used  to  mark  apparent  unidirectional  segments.  The 
squared  multiple  coherence  analysis  at  WMO  indicated  little 
if  any  unidirectional  noise,  and  no  dominant  direction  is 
apparent  on  the  filtered  seismogram.  The  data  used  is  about 
the  same  as  that  for  Figure  11. 

WMO  strain  seismogram  rotation 

In  order  to  rotate  the  strain  seismograms  it  is  necessary 
to  use  plane  strain  theory  which  implies  that  two  of  the  princi¬ 
pal  strain  axes  are  in  the  horizontal  plane.  It  is  also 
necessary  to  use  a  third  strain  seismometer  to  determine  the 
shear  component  of  strain.  Jaeger  (1962,  p.  42)  writes  the 
rotated  components  of  strain  relative  to  the  original  components 
in  the  form 


2q 

=  ex  cos  9 

+  Y  sin9 
'  xy 

cos  6  + 

2 

e^sin  9 

2 

=  e  sin  9 

-  y  sin9 

cos8  + 

2 

e  cos  9 

X 

'  xy 

y 
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In  which  ,  £y  are  the  linear  strain  components  and  y 
is  the  shear  component  of  the  strain  tensor,  and  8  is  £he 
angle  through  which  the  ..train  coordinates  are  rotated.  Desig¬ 
nating  the  north  component  of  strain  by  SN,  the  east  by  SE 
and  the  northwest  by  SNW,  the  above  equations  becomes 

^  2 

ex  =  SEcos  6  +  ( SE+SN- 2SNW)sin0cos0  +  SNsin2© 

^  2 

Ey  =  SEsin  6  -  (SE+SN- 2SNW)sin6cos0  +  SNcos20 

The  45  degrees  azimuthal  difference  between  north  and  north¬ 
west  gives  this  equation  a  particularly  simple  form. 

Figures  23  through  26  are  the  result  of  applying  the 
above  formula  to  three  of  the  horizontal  strain  components  at 
WMO.  Apparent  unidirectional  segments  of  data  have  been 
enclosed  in  rectangles.  The  analysis  could  not  be  repeated 
at  HNME  because  of  the  absence  of  horizontal  strain  seismometers 


Squared  coherence  of  short  samples  at  WMO 

Figures  27  and  28  are  plots  of  the  squared  coherence 
between  the  vertical  strain  and  vertical  pendulum  at  WMO. 

The  samples  occur  in  two  groups.  Analysis  of  the  squared 
coherence  values  and  the  sample  intervals,  show  that  the 
relatively  high  multiple  coherence  of  one  sample  is  caused 
by  the  coherence  at  the  data  between  points  813  and  913. 

Figure  29  is  a  copy  of  the  original  seismogram  on  which  the 
coherent  segment  has  been  marked.  Examination  of  the  rotated 
horizontal  pendulums  for  this  interval  (Figure  20)  shows  that, 
the  azimuthal  direction  of  the  coherent  noise  segment  is 
90  degrees  (or  270  degrees),  along  the  axis  of  the  Wichita 


* 
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Mountain  thrust.  The  rotated  strain  seismograms  do  not  exhibit 
a  well-defined  azimuthal  direction.  There  exists  a  pulse  from 
1Q5  degrees  Cor  285  degrees)  which  may  be  the  same  noise  pulse 
as  seen  on  the  horizontal  pendulum  instruments.  Figure  29 
shows  the  90  degrees  phase  shift  expected  for  Rayleigh  waves. 

Noise  at  HNME  at  2.6  Hz 

At  HNME  for  all  samples  there  exists  a  peak  in  the  noise 
spectra  at  2.6  Hz.  The  noise  has  a  sinusoidal  appearance  on 
most  seismograms.  A  number  of  lines  of  evidence  described 
below  show  that  it  is  P-wave  noise.  The  mechanism  of  its 
generation  is  left  for  later  study. 

Figure  30  is  a  detailed  plot  of  the  relative  power  on  the 
vertical  strain  and  vertical  pendulum  instruments  at  HNME.  The 
higher  peak  of  the  pendulum  relative  to  the  strain  is  indica¬ 
tive  of  P-wave  noise,  since  the  vertical  strain  is  less  responsive 
to  P-waves  than  the  vertical  pendulum  (Romney,  1964).  This  same 
relative  amplitude  phenomena  is  also  characteristic  of  the  2.0  Hz 
noise  at  WMO. 

Figure  31  is  the  first  motion  of  a  well-defined  earthquake 
event  at  HNME,  which  shows  that  PZ  and  SZ  are  in  phase  for 
P-waves.  Figure  32  is  the  phase  plots  from  four  samples  which 
show  a  phase  difference  of  360  degrees  (or  0  degrees)  at  2.6  Hz. 

The  squared  multiple  coherence  of  points  500-2547  HNME 
(Seis.  19028)  has  a  value  of  0.92  at  2.6  Hz  (Figure  11).  This 
means  that  for  this  sample  the  noise  is  nearly  unidirectional. 
Other  samples  have  lower  squared  multiple  coherence  indicating 
that  the  unidirectionality  of  this  sample  is  not  typical  and 
therefore  should  not  be  ascribed  to  cultural  noise  originating 


at  a  particular  azimuth.  Figure  33  are  the  rotated  seismograms 
corresponding  to  the  above  described  squared  multiple  coherence. 
The  azimuth  is  well-defined  at  168  degrees  (or  348  degrees).  The 
cosine  azimuthal  response  of  the  horizontal  instruments  makes 
thi*  value  only  approximate  (say  about  15°).  The  seismograms 
of  Figure  33  have  been  band-pass  filtered  1.3  to  5.2  Hz,  24 
db/octave  phaseless. 

34  consists  of  samples  from  the  168  degree  azimuth 
plot  of  Figure  33  compared  with  the  vertical  pendulum  filtered 
1.3  to  5.2  Hz.  These  samples  show  the  general  0  degree  phase 
differences  expected  from  the  above  evidence. 

Noise  at  HNME  0.5  to  1.2  Hz 


At  HNME  the  squared  coherence  vertical  strain  to  vertical 
pendulum  is  commonly  about  0.9  between  0.15  and  0.5  Hz,  which 
means  that  there  exists  a  well-defined  linear  transformation 
that  can  be  used  to  cancel  noise  in  this  frequency  bc'nd.  Between 
0.5  and  1.5  Hz  the  squared  coherence  decreases  to  values  less 
than  0.2  which  is  below  the  95  percent  confidence  limit  that  the 
linear  transfer  function  in  the  frequency  domain  is  not  zero 
(Figure  8).  Nonetheless  the  estimated  phase  angle  remains  relatively 
constant  at  about  270  degrees  from  0.15  to  1.5  Hz.  (Figure  32).  This 
may  reflect  the  existence  of  several  different  independent  Rayleigh 
modes,  each  with  a  different  transfer  function  amplitude  but  identi¬ 
cal  270°  phase  shift. 

Rotation  of  the  horizontal  pendulum  instrument  exhibits  a  few 
short  segments  of  noise  even  between  0.5  and  1.5  Hz  which  appear 
to  be  unidirectional.  In  most  cases  the  noise  is  not  thus 
resolvable,  meaning  that  the  noise  is  arriving  simultaneously 
at  the  site  from  two  or  more  azimuths.  Figures  35,  36,  and  37 
are  examples  of  apparent  unidirectional  noise  from  Seismogram 
19028  between  points  500  and  2500.  These  data  are  also  plotted 


m  F igure  16  while  figures  15  and  17  give  the  corresponding 
plots  for  the  frequency  hands  0.1  to  Q.5  and  1.2  to  2.4  Hz. 

The  criteria  for  a  single  direction  noise  is  low  noise  level 
on  the  instrument  at  90  degrees  to  the  apparent  source  azimuth. 

The  portion  of  noise  analyzed  has  been  band-pass  filtered 
(0.5  to  1.2  Hz,  24  db/octave,  phaseless). 

Figure  35,  the  shaded  pulse,  shows  a  direction  whose 
azimuth  is  168  degrees  (or  348  degrees).  The  vertical  pen¬ 
dulum  is  90  degrees  out  of  phase,  and  the  vertical  strain  is 
180  degrees  out  of  phase  with  the  rotated  horizontal  seismogram 
as  would  be  expected  for  Rayleigh  waves . 

Figure  36,  the  shaded  pulse,  illustrates  a  unidirectional 
noise  pulse  that  is  apparently  not  a  Rayleigh  wave.  There  should 
be  a  0°  or  180°  phase  difference  between  vertical  strain  and 
horizontal  pendulum  for  Rayleigh  waves.  This  noise  pulse  has 
a  90  degree  phase  difference.  The  fact  that  the  vertical 
pendulum  is  180  degrees  out  of  phase  with  the  rotated  h  izontal 
pendulum  indicates  a  P  or  SV— wave.  The  90  degree  phase  differ¬ 
ence  between  the  vertical  strain  and  vertical  pendulum  also  indi¬ 
cates  a  P  or  SV  wave.  If  it  is  an  SV-wave  it  may  be  shown  that 
its  angle  of  incidence  is  less  than  critical. 

Figure  37  illustrates  two  pulses  on  the  horiztonal  instru¬ 
ments  which  seem  to  merge  into  a  single  pulse  on  the  vertical 
instruments.  The  first  portion  (3  degree  azimuth)  has  the 
characteristics  of  a  Rayleigh  wave.  The  93  degree  azimuth 
pulse  is  in  phase  on  all  three  instruments.  Preliminary  identi¬ 
fication  is  that  this  is  an  SV-wave  whose  angle  of  incidence  is 
greater  than  the  critical  angle. 

This  brief  examination  of  the  noise  field  in  the  0.5  to  1.2  Hz 
band  suggests  that  the  decreasing  squared  coherence  above  about 
0.5  Hz  is  due  to  a  mixture  of  wave  types. 
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CONCLUSIONS  AND  CURRENT  WORK 


At  HNME  the  evidence  is  clear  that  short  period  noise, 

0.15  to  0.60  Hz,  can  be  cancelled  by  a  linear  combination  of 
the  vertical  strain  and  vertical  pendulum.  Over  80  percent  of 
this  noise  can  be  shown  to  be  unidirectional.  Between  0.6  and 
2.0  Hz,  the  coherence  is  low,  indicating  that  no  linear  combi¬ 
nation  exists  which  will  cancel  the  noise.  This  noise  apparently 
comes  from  many  directions.  Despite  this,  the  relative  phase 
between  the  two  instruments  is  stable  at  90  degrees  over  the 
entire  frequency  band  0.1  to  2.0  Hz,  indicating  perhaps  that 
some  small  component  of  the  noise  between  0.6  and  2.0  Hz  is 
Rayleigh  waves  which  may  be  cancelled  out.  An  alternative 
hypothesis  is  that  of  energy  in  several  different  Rayleigh  modes.  ' 
At  2.6  Hz  there  is  a  well  defined  band  of  P-wave  noise. 

Further  studies  at  HNME  are  concentrating  on  analysis  of 
the  several  phases  generated  by  an  earthquake  near  Morocco,  with 
a  view  to  seeing  how  each  of  them  is  detected  on  the  complete 
array  of  instruments. 

WMO  has  no  frequency  band  in  which  a  well-defined  transfer 
function  exists  or  in  which  the  noise  is  unidirectional.  Further 
research  will  be  pointed  toward  giving  a  physical  explanation 
of  why  this  is  in  order  to  establish  criteria  for  the  applicability 
of  the  strain  system  to  a  particular  site.  We  plan  to  do  this  by 
examining  earthquake  phases  which  arrive  at  WMO  from  different 
c.'imuths  traversing  nearby  geologic  structures  of  greater  or 
ltsser  complexity. 
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Figure  13.  HNME  -  Squared  Multiple  Coherence. 
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Figure  34.  HNME  -  Ph«se  Comperlson  -  2.6  Hr. 


Figure  35.  HNME  -  Phase  Comparison  (0.5-1. 2  Hz). 


APPENDIX  I 


Data  Used  in  Preparation  of  Report 

Seis.  18403  WMO  16,000  points  beginning  at  0638  27  January 
1969  -  includes  earthquake  in  Soviet  Arctic  East  of 
Severnaya  Zemlya. 

Seis.  18736  HNME  16,000  points  beginning  at  0520  28  February 
1969  noise. 

Seis.  18767  FOX  16,000  points  beginning  0452  22  November 

1968  noise. 

Seis.  18774  HNME  16,000  points  beginning  0425  28  February 

1969  includes  N.  Atlantic  earthquake. 

Seis.  19028  HNME  16,000  points  beginning  0600  15  March 
1969  noise. 

Seis.  19029  HNME  16,000  points  beginning  1337  15  March 
1969  includes  Aleutian  earthquake. 
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